Electron correlations and single-particle physics in the Integer Quantum Hall Effect 
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The compressibility of a two-dimensional electron system with spin in a spatially correlated ran- 
dom potential and a quantizing magnetic field is investigated. Electron-electron interaction is treated 
with the Hartree-Fock method. Numerical results for the influences of interaction and disorder on 
the compressibility as a function of the particle density and the strength of the magnetic field are 
presented. Localization-delocalization transitions associated with highly compressible region in the 
energy spectrum are found at half-integer filling factors. Coulomb blockade effects are found near 
integer fillings in the regions of low compressibility. Results are compared with recent experiments. 



The integer quantization of the Hall conductance of 
a two-dimensional electron system (2DEG) in a strong 
magnetic field can be understood in terms of quan- 
tum phase transitions near the centres of the Lan- 
dau bands associated with disorder-induced localization- 
delocalization transitions of single-electron states in a 
one parameter scaling model. Neglecting interaction, 
the localization length has been found diverging, £ oc 
\E — Ec\~ h ', with the critical energy E c . The universal 
value of the critical exponent, v = 2.34 ± 0.04 [HQ- is 
widely accepted. Peaks in the magneto-conductance are 
associated with E c . The localized states in the band tails 
are associated with zero conductance at zero temperature 

a. 

Despite consistency with several transport experiments 
0, the validity of the one-parameter scaling model has 
been controversely discussed, in particular the impact 
of interaction on the frequency-dependent scaling of the 
conductivity and the tunneling density of states jy, 0, 

In recent experiments, mesoscopic conductance fluctu- 
ations found in silicon metal-oxidc-scmiconductor field- 
effect transistors in dc-transport show regular patterns 
which have been interpreted as charging effects |9(. Pat- 
terns associated with Coulomb blockade in localized 
states have also been found in measurements of the shift 
d\ij dn of the chemical potential \i with a scanning single- 
electron transistor probe when changing the particle den- 
sity n |10|. The latter results have been interpreted in 
a model in which the quantum Hall transition appears 
as a result of the strong and complete screening of the 
disorder potential near half-integer filling factors. The 
absence of screening in the incompressible regions of the 
energy spectrum leads to localized states that account for 
the observed charging effects. In this model, the phase 
transition has been interpreted as a percolation transi- 
tion between incompressible and compressible regions at 
certain concentrations of localized charge islands. 

We report results of an extensive unrestricted Hartree- 
Fock (HF) study of the 2DEG with spin in the presence of 
long-range correlated disorder and a perpendicular mag- 
netic field that can contribute towards more detailed un- 
derstanding of the experiments. We have studied the 
change d/i/dn of the chemical potential /x with the par- 



ticle density n = N/L 2 of N particles in a square sys- 
tem of linear size L as a function of n and the magnetic 
field strength B. This quantity is proportional to the in- 
verse compressibility k~ 1 oc d^i/dn. We find interaction- 
induced enhancement of the g-factor and strong evidence 
for charging effects in regions near integer filling factors. 
The shapes of the HF quasiparticle wave functions in- 
dicate no significant change in the localization behavior 
at the Fermi level compared to the non-interacting limit. 
This is due to quantum corrections which modify the per- 
colation mechanism [ill. O ur results are consistent with 
the recent experiment |10|. However, we find that the 
interaction does not destroy the critical behaviour at the 
quantum Hall phase transition. 

The Hamiltonian of the 2DEG in GaAs is H° + V c 
with = (p - eA) 2 /(2m*) + sg^ B B/2 + V dis (r) with 
the vector potential A ~ (0,Bx,0), flux density B 
and spin s = ±1. m* = 0.067m o is the electron 
mass, fiB is the Bohr magneton and g = —0.44 the 
electron ^-factor. The impurity potential is Vd; s (r) = 
E^i^/^WKr-rOVd 2 ] with JVj the number of 
scatterers at random positions with random strengths 
Vi, —Vq < Vi < Vq. The range d of the impurity po- 
tential is the spatial correlation length of the random- 
ness, d — corresponds to uncorrelated disorder, d > Ib 
(Ib = v h/m*uj c magnetic length, lo c — eB/m cyclotron 
frequency) yields a slowly varying potential which is be- 
lieved to be adeqate for high mobility samples. The disor- 
der introduces the energy scale V = (N\Vq /IgL 2 ) 1 / 2 (L 2 
area of the 2DEG). The Coulomb interaction Vc(r — r') = 
e 2 / (4-7ree|r— r'|) introduces an energy scale 7 = e 2 /47ree^ 
(e elementary charge, dielectric constant e = 12.4 ). Peri- 
odic boundary conditions are assumed 0, 0] . Neglect- 
ing disorder and interaction, the Schrodingcr equation 
yields the Landau wavefunctions \mX) (X = kjl 2 B guid- 
ing center coordinate, kj = 2nj / L wavenumber) that are 
used for the construction of the HF basis. 

The HF equation is 
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E as c o 



where (mX\as) = C^ x = C" s are the expansion coef- 
ficients of the HF states \as) and E as the energy eigen- 



2 



values. The Fock matrix F^ Xm , x , = F£ = Hfa + 
E E PabMijab - pl b M iabj has to be determined selfcon- 

a b 

sistently. It contains the interaction matrix elements 

Mijab = £ E V(q)(i\e^\]}{a\e-^\b) and the density 

q*,q v 

matrix p«=E.Pfi= E C^*Cf. 

a(occ) 

We have performed selfconsistent calculations for B in 
the range of 1 . . . 6T and n — (2 . . . 200)/L 2 in a square 
of length L = 30a for a = 10 nm. The sample size is 
comparable to the sample size studied in Ref. E This 
yields electron densities n = (0.22 ... 22) 10 10 cm" 2 . We 
used 250 impurities with d = 2a and Vq = 2 meV to 
model a high mobility sample. 

To study Coulomb interactions, it is useful to control 
its strength relative to the other energy scales. We set 
7 = C70; c = corresponds to the noninteracting system, 
c = 1 to the non-screened Coulomb interaction. For each 
combination of parameters, the self-consistent field and 
thus total energy and chemical potential are determined. 
Thus, the self-consistent field can relax with respect to 
changes in parameters and ground state and the quasi- 
particle wavefunctions are optimized for given N and B. 
This is sometimes referred to as the "dclta-SCF- method" 
[lij : The difference between total energies is calculated 
successfully in HF approximation, because the response 
of the iV-particle system to the addition of another elec- 
tron or hole is contained in the calculation of the N ± 1- 
particle many-body wavefunction. 

The number of Landau bands per spin included in the 
calculation was chosen such that if initially the electron 
with the highest energy at the highest particle num- 
ber (N max = 200) is in Landau band m max , m max + 1 
Landau bands per spin are used in order to provide 
enough states for the self-consistent field calculations; 

= L 2 / (2nl 2 B ) = BL 2 e/h is the number of flux quanta 
in the system and equals the number of states per spin in 
a Landau band. For example, at B = 1 T and = 21, 
we have used 11 Landau bands per spin, in total 462 
states to host 2 to 200 electrons. In the results presented 
below, we focus on particle numbers 2 < N < 117. 

For energies in the band tails, the wave functions are 
localized near equipotential lines, at least for strong mag- 
netic field. In the center of the band, the wave functions 
are delocalized (Fig. [I] inset). 

The scaling of the participation number P at the Fermi 
energy E, P^ 1 = j d 2 r|V> 4 (r)| 4 oc £~ 2 (E) as a func- 
tion of the filling factor (Fig. ^) has been investigated 
[lfij . With i> ~ 2.3 rcsonable collapse to a single curve 
has been achieved of the data for different system sizes 
consistent with the critical behavior of the localization 
length without interaction. The curve resembles the scal- 
ing of HF wave functions obtained at fixed filling factor 
using occupied and empty HF orbitals 0]. Within our 
method the orbitals are separately determined for ev- 
ery combination of electron density and magnetic field, 
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FIG. 1: Scaling of the participation ratio II = L 2 ~ D ^ 2 ' P as a 
function of the filling factor (y — ^ C )L 1// ", disorder-averaged 
for each system length L flEf . Parameters v = 2.3 (critical 
exponent) and D(2) = 1.6 (correlation dimensions of wave- 
functions |2"H0 Lead to reasonable data collapse. Inset: De- 
pendence of P v (Ef) on v for spin-up (solid) and spin-down 
(dashed) for L = 30.1a, disorder-averaged. The g- factor en- 
hancement is so large, that a subsequent filling of the spin- 
split Landau level is observed. Color inset: Typical probabil- 
ity density [red (bright) high, blue (dark) low value] of a HF 
state near the band center. 



together with the self-consistent potential. The consis- 
tency within the errors of the data with the scaling hy- 
pothesis indicates that the critical behavior of the HF 
states is unaffected by the change of the self-consistent 
field induced by changing the filling factor. However, 
this must be confirmed by more precise scaling studies 
for larger systems Q. The HF energy for N particles 
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chemical potential /j, = E^~ l ' B — E^ B and 
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This definition reflects a global property, in contrast to 
the local compressibility reported in Ref. However, 
the total energy and its derivative depend on the density 
matrices p''* formed with the HF orbitals. The latter 
are obtained for a finite system and can be considered to 
reflect local properties such as a specific disorder poten- 
tial and electrostatic and exchange interaction with the 
surrounding electrons. Thus, the inverse compressibility 
calculated here can be expected to reflect the features 
observed in the measurements done with a tunnel tip. 

Figure[3]shows dfx/dn for various particle numbers and 
magnetic fields without and with interaction. Without 
interaction (Fig. |2t), states are compressible almost 
everywhere except near even integer filling, v = N /N^. 
This is expected since every particle that enters finds 
many states at energies close to the Fermi level. The 
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FIG. 2: d^/dn (in units of meV/L 2 ) in the (n, _B)-plane 
for interacting electrons (a) 7/70 = 0, (b) 7/70 = 0.1, (c) 
7/70 = 0.5. Dark (blue) regions high, bright (red) regions 
low compressibility. Solid white lines are guide to the eyes. 



Fermi energy remains almost constant with respect to 
the particle number, as long as the actual Landau band 
is not completely filled. Near even integer filling, the 
Fermi level jumps to the next Landau level due to very 
small level density. In the presence of Zeeman splitting, 
there are N$ states per spin. Since the Zeeman energy in 



GaAs is very small as compared to other energy scales, 
the levels overlap strongly and are equally occupied. In- 
compressible lines are then only obtained for even fill- 
ings. With interaction, the Zeeman splitting is large due 
to exchange enhancement of the g- factor |l6j . Spin-up 
and spin-down levels are well separated, resulting in ad- 
ditional incompressible lines at odd fillings. This effect is 
exaggerated in the unrestricted HF approximation . 

Figures[2b>c show compressibility patterns for different 
interaction strengths. The compressibility of the inter- 
acting electrons shows several regular structures. Hor- 
izontal lines of constant compressibility parallel to the 
B-axis appear below v = 1/2, enclosed by lines of low or 
even negative compressibility. At v w 1/2, there is a re- 
gion of high compressibility. For filling v w 1 , lines of low 
compressibility parallel to n = vtlb = j /2irl B = jeB/h 
are observed. The range of electron density where this 
happens is independent of B. The number of the strongly 
localized states must therefore be independent of B. This 
has been ascribed before to Coulomb blockade in stro ngly 
localized states associated with deep potential wells [l(| . 
When the potential landscape is completely screened the 
addition of a further electron is possible. This can even 
result in a negative compressibility, which in this case is 
not related to a thermodynamic instability, but to the 
fact that a positive impurity can be overcompensated by 
an entering electron. This in turn causes a depletion of 
electronic charge afterwards in this region [18j. These 
charging effects occur also in the higher Landau levels, 
although less prominent, because the localization length 
of the electrons increases with the Landau level index, 
thus making the system more sensitive to compression. 
We emphasize that the behavior of the HF total energy 
with n and B provides the correct charging pattern in 
the strongly localized region. 

In addition to these Coulomb interaction-dominated 
features, we observe highly compressible regions around 
half integer fillings v = j/2, Ub = jeB/2h. These cor- 
respond to the centers of the Landau bands where the 
disorder is considerably screened. The width of these 
regions, An ex t, is roughly constant as a function of B. 
The number of the effectively extended states (diame- 
ters larger than L) in a Landau band must be almost 
independent of B although the total number of single- 
electron states per Landau band increases linearly with 
B. This can be understood in the one-band approxima- 
tion, where the single particle density of states D scales 
as D(E/T) = (N < j,/T)f{EB/Y) [H|. The energy interval 
AE = \E — E c \ in which the localization length exceeds 
the system size is defined by £(E) — £o\E — Eo\~' y > L. 
Thus, AE = (L/Zo)- 1 / oc B- 1 /", since £ oc T~ 2 oc 
B~ l , and Sn cxt w D(0)AE oc B 1 / 2 " 1 / 5 = B 065 with 
v k, 2.3. 

In contrast to previous assertions [lflj . we find be- 
tween the low-compressibility regions of Coulomb block- 
ade in the strongly localized states and the high com- 
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prcssibility regions of derealization, that there arc large 
regions of intermediate statistically fluctuating compress- 
ibility. These correspond to localized states that cover 
larger spatial regions with randomly fluctuating areas. 
The charging energies of these states, if applicable at all, 
should be much smaller than in the regimes of strongly 
localized states, and also strongly fluctuating. As a con- 
sequence, one would not expect regular compressibility 
patterns in these intermediate regions, and this is what 
is observed in Fig. In these regions, the localization 
properties arc determined by the competition of tunnel- 
ing between, and destructive interference along the perco- 
lating cquipotcntial lines, and it is this competition that 
is responsible for the critical behavior [Jjj. The regimes 
of strong localization and extended states are clearly ob- 
served, separated by regions of intermediate states. 

This could be affected by the finite system size L, if the 
localization length of the intermediate states is compara- 
ble to L. Here, it is likely that the fluctuations vanish in 
the thermodynamic limit, and only incompressible local- 
ized stripes remain, consistent with experiment. Prelimi- 
nary results for larger systems show that the fluctuations 
seem to remain present. A systematic study of the size 
dependence will be reported elsewhere |2?i| . 

In conclusion, we have investigated the density depen- 
dence of the chemical potential as a function of electron 
density and magnetic field for a quantum Hall system. 
We have shown that electron interactions, treated in HF 
approximation, but with the possibility for the ground 
state to respond to changes in magnetic field or electron 
density, modify the compressibility pattern. The appear- 
ance of regular structures can be interpreted as charging 
of localized states. This is consistent with recent ex- 
periments and suggests that interactions are important 
for the understanding of the integer quantum Hall effect, 
especially in the plateau regions. However, the results 
reported here are not in contradiction to the conjecture 
that the critical behavior of the metal-insulator transi- 
tion is unaffected by interactions and microscopic details 
of the disorder potential since the scaling of the partici- 
pation number at the Fermi energy is found to be consis- 
tent with the scaling hypothesis, although for each elec- 
tron density the effective potential changes as a result of 
charge rearrangement. 

The question remains open if correlations beyond the 
HF approximation can affect the compressibility pattern. 
Calculations regarding correlation effects have not been 
done. However, the HF results reproduce the charging 
effects in the regimes of localized states well and support 
the assumption that the critical behaviour of the inte- 



ger quantum Hall transition can be understood within a 
single particle model. 
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